Method for Initialization and Stepsize Control of 
Time-Domain Equalizer in Multi-carrier Communication System 

Background of the Invention 

5 

(a). Field of the Invention 

The present invention relates in general to a multi-carrier 
communication system, and more particularly to a method for initiation and 
stepsize control of a time-domain equalizer in a multi-carrier 
10 communication system. 

(b^i. Description of the Prior Arts 

Muhi-carrier modulation is widely used in communication systems 
nowadays. Fig.l is a block diagram of a communication system 100 using 
multi-carrier modulation. The communication system 100 employs a set of 

15 N-point Inverse Fast Fourier Transform (IFFT) 102 in the transmitter and 
the Fast Fourier Transform (FFT) 111 in the receiver to transceive data. A 
channel 106 is divided into N subchannels for transmitting data. The signal 
transmitted by one subchannel is orthogonal to those transmitted by any 
other subchannels. Thus, the signal transmitted by the subchannels would 

20 not interfere each other, and inter-channel interference (ICI) can be avoided. 

The set of N-point data outputted from IFFT 102 is called a symbol. 
Since the channel impulse response (CIR) is not ideal, that is, the 
amplitude/frequency response of the chaimel are not constant across the all 
used subchannels, then the received signals will different from the 

25 transmitted signal, and the signals presented to the QAM decoder 113 will 
different from the signals outputted from the QAM encoder 101. If the 
distortion is severe then the data transmitted by one subchannel will affect 
both the data transmitted by other subchannels (ICI) and the data transmitted 
by that subchannel in the previous and subsequent symbol periods 

30 (inter-symbol interference, ISI). In order to avoid inter-symbol interference 
(ISI) and ICI, a "cyclic prefix"(CP) is added to each symbol, i.e. the last u 
points of each symbol are copied and added in the front of the symbol, as 
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shown in Fig.2. Therefore, each symbol outputted from adding cyclic prefix 
circuit 103 includes (N+v) points. The adding cyclic prefix circuit 103 in the 
transmitter and the removing cyclic prefix circuit 110 in the receiver of Fig. 1 
are used to add and remove cyclic prefixes respectively. 

5 In the communication system 100 of Fig.l, if the valid length of 

channel impulse response (CIR» denoted by h[n]) is shorter than the length 
of cyclic prefix, then a symbol, after being transmitted in the channel 106 
and received by the receiver (i.e. the convolution of the symbol and the CIR, 
h[n]), would not affect the data of the subsequent symbol received in the 

10 subsequent symbol period. However, if the length of CIR is larger than that 
of the cyclic prefix, then ICI and ISI will occur. Under this circumstance, a 
time-domain equalizer (TEQ) with an impulse response w[n] ,as shown in 
the block 108 of Fig.l, is necessary for the receiver of the system 100. The 
TEQ 108 is used to modify the CIR of the communication system 100, such 

15 that valid length of the modified CIR (called target impulse response and 
denoted by b[n]), i.e. the convolution of the CDR. h[n] and the TEQ impulse 
response w[n], is shorter than that of the cyclic prefix, thereby preventing 
the received data firom ICI and ISI. 

Since the CIR is different with various transmission channels, the TEQ 
20 impulse response needs to be adjusted accordingly. Many adaptive TEQ 
algorithms are developed in succession. Since these adaptive TEQ 
algorithms are sensitive to initial values of TEQ impulse response W[n] and 
target impulse response B[n] in frequency-domain (In the following 
description, time-domain variables will be referred to by lower-case letters 
25 and frequency-domain variables will be referred to by capital letters), the 
adapting result would be unreliable if the initial values are not properly 
determined dimng the TEQ Initialization process. The adapting result will 
fail into the local maxima point. 

U.S. Pat. No. 5,285,474 and 6,396,886 disclosed the conventional 
30 adaptive TEQ algorithms. However, both of these two patents failed to 
disclose the method to determine the initial values of frequency-domain 
TEQ impulse response W[n] and target impulse response B[n]. In addition, 
the coefficient (called stepsize coefficient in this specification) used to 
stepwise adjust W[n] using frequency-domain Least Mean Square (LMS) 
35 will greatly affect the performance of the adaptive TEQ algorithms. If the 



stepsize coefficient is set too small, the converging speed will be too slow; if 
too large, then a diverging result will occur frequently. Both cases would 
degrade the system performance seriously. 

5 Sammarv ftf the Invention 

In view of the above-description, it is one of the objects of the present 
invention to provide a method for initializing and stepsize controlling a 
time-domain equalizer (TEQ) for use in a muhi-carrier communication 
10 system to upgrade the performance of the TEQ performing the adaptive 
TEQ algorithms. 

According to the object of the present invention disclosed above, one 
of the embodiment of the present invention discloses a method for 
initializing a time-domain equalizer (TEQ) comprised in a receiver of a 

15 muhi-carrier commimication system, the method comprising: estimating a 
channel impulse response (CIR) h[n] according to a received symbol, 
wherein the received symbol includes a cyclic prefix with v points and a 
data portion with N points, wherein i=0~N-l; selecting one of a plurality of 
groups according to the total energy of the groups, wherein each group 

20 includes consecutive v points of the received symbol; modifying the channel 
impulse response (CIR) h[n] according to the selected group to generate a 
modified channel impulse response h'[n]; generating a target impulse 
response b[n] according to the modified channel impulse response h'[n] and 
a window mask m[n], wherein i=0~N-l; transforming the channel impulse 

25 response (CIR) h[n] and the target impulse response b[n] to a 
frequency-domain to generate a frequency-domain channel impulse 
response H(i) and a frequency-domain target impulse response B(i) 
respectively, wherein i=0~N-l; and generating a frequency-domain impulse 
response W(i) of the time-domain equalizer according to the 

30 frequency-domain chaimel impulse response H(i) and the frequency-domain 
target impulse response B(i). The frequency-domain target impulse response 
B(i) and the frequency-domain TEQ impulse response W(i) are for 
initializing the TEQ. 
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Another embodiment of the present invention discloses a method for 
initialising a time-domain equalizer (TEQ) comprised in a receiver of a 
multi-carrier communication system, the method comprising: estimating a 
channel impulse response (CIR) h[n] according to a received symbol, 

5 wherein the received symbol includes a cyclic prefix with v points and a 
data portion with N points, wherein i=0-N-l; selecting one of a plurality of 
groups according to the total energy of the groups, wherein each group 
includes consecutive v-lw points of the received symbol, wherein Iw is the 
length of the TEQ impulse response; modifying the channel impulse 

10 response (CIR) h[n] according to the selected group to generate a modified 
channel impulse response h'[n]; determining a frequency-domain impulse 
response W(i) according to a frequency-domain modified channel impulse 
response H^i), wherein i^O-'N-l; and determining a frequency-domain 
target impulse response B(i) according to the frequency-domain impulse 

15 response W(i) and a frequency-domain CIR H(i), wherein i=O^N-l. The 
frequency-domain target impulse response B(i) and the frequency-domain 
TEQ impulse response W(i) are for initializing the TEQ. 

Still another embodiment of the present invention discloses a method 
for adapting a time-domain equalizer (TEQ) comprised in a receiver of a 

20 multi-carrier communication system, the method comprising: determining a 
frequency-domain TEQ impulse response W^(i) and a frequency-domain 
target impulse response Bk(i) for initializing the TEQ, wherein i^O-^N-l; 
generating a modified TEQ impulse response wi^w(i) and a modified target 
impulse response hj^y^(i) according to a time-domain TEQ impulse response 

25 wk(i) and a time-domain target impulse response bk(i); determining a error 
term Ek(i) according to the modified TEQ impulse response wic,w(i)j a 
modified target impulse response bicw(i)> and a frequency-domain channel 
impulse response (CIR) H(i) ; adjusting the frequency-domain TEQ impulse 
response Wk(i) to generate a adjusted frequency-domain TEQ impulse 

30 response Wk+i(i) through performing a least mean square (LMS) operation 
according to the error term Ek(i), a frequency-domain modified TEQ 
impulse response Wic,v»^(i), a frequency-domain channel impulse response 
(CIR) H(i), and a stepsize coefficient ji, wherein the stepsize coefficient \i in 
a time-varying coefficient; generating a modified adjusted 

35 frequency-domain TEQ impulse response Wic-ri,w(i) according to the adjusted 
frequency-domain TEQ impulse response Wit+i(i); and adjusting the 



frequency-domain target impulse response Bk(i) to generate a adjusted 
frequency-domain target impulse response Bk+i(i) according to the modified 
adjusted frequency-domain TEQ impulse response Wjc+i,w(i) and the 
frequency-domain channel impulse response (CIR) H(i). 

5 These and other objectives of the present invention will no doubt 

become obvious to those of ordinary skill in the art after reading the 
following detailed description of the embodiments that is illustrated in the 
various figures and drawings. 

10 Brief Description of the Drawings 

Fig.l is a block diagram of a communication system using 
multi-carrier modulation. 

Fig,2 is a diagram showing the formation of a cyclic prefix. 

15 Fig.3 is a flow chart of the method for determining the initial values of 

frequency-domain TEQ impulse response W[n] and target impulse response 
B[n] according to the first embodiment of the present invention. 

Fig.4A and 4B are diagrams showing two examples of the time-domain 
window mask. ^ 

20 Fig-5 is a flow chart of the method for determining the initial values of 

frequency-domain TEQ impulse response W[n] and target impulse response 
B[n] according to the second embodiment of the present invention. 

Fig.6 is a flow chart of a typical adaptive TEQ algorithm. 

Fig.7. is a diagram showing an example of the fimction f(k) used to 
25 determine the stepsize coefficient function according to the embodiment of 
the present invention. 

Detailed Description of the Embodiments 
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Fig. 3 is a flow chart of the method for determining the ini tial values of 
frequency-domain TEQ impulse response W[n] and target impulse response 
B[n] according to the first embodiment of the present invention. As 
mentioned above, the TEQ 108 is used to shorten channel impulse response 
5 (CIR) of the multi-carrier communication system 100. The CIR is 
represented by a N-point sequence h[n], n=0-N-l. The communication 
system 100 encodes data as symbols for transceiving, wherein each symbol 
comprises N sampling points and a cyclic prefix with v sampling points. 

As shown in Fig.3, the first embodiment for determining the initial values of 
10 frequency-domain TEQ impulse response W[n] and target unpulse response 
B[n], which is called the Time-domain Window Mask method in the 
specification, includes steps of: 

30 producing an estimated CIR h[n] according to the received 
1 5 symbols which are known by both of the receiver and the transmitter; 

31 calculating the total energy of each group of consecutive v 

y 

points and selecting the group of consecutive v points with the 
maximum energy; 

32 setting all the remaining points outside the selected group to zero 
20 to generate a modified estimated CIR h[n] ; 

33 multiplying the modified estimated CIR h[n] with a time-domain 
window mask m[n] to generate an estimated target impulse response 

b[n], wherein b[n]=h[n] • m[n], n=0~N-l; 

34 transforming h[n] and b[n] into the frequency domain to generate 
25 H(i) and B(i) respectively, i-0~N-l, wherein B(i) is determined to be 

the initial value of the frequency-domain target impulse response; and 

35 determining W(i) by dividing B(i) by H(i), i=0~N-l, wherein W(i) 
is determined to be the initial value of the frequency-domain TEQ 
impulse response. 
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The procedure of step 31 that the v consecutive points with the 
highest total energy are found by a cyclic search through N points of h[n], 
that is, X|/2[«]|^ , k=0'"N"l, is called the Locate Maximum Energy 

algorithm. The number of consecutive pomts of h[n] to perform the Locate 
5 Maximum Energy algorithm is determined such that the valid length of the 
estimated target impulse response b[n] is shorter than that of the cyclic 

prefix (having v points). In this embodiment, it is v consecutive points 
of h[n] with the highest total energy to be found through performing the 
Locate Maximum Energy algorithm. However, the value v is not the 
1 0 limitation of the present invention. 

In the step 33, if the points of h[k] - h[k+ v -1] are selected, the points 

of m[k] - m[k-»- v -1] of the time-domain window mask m[n] are set to be 

the value between zero and one and the other points of m[n] are set to be 
zero to further modify the modified estimated CIR h[n3. 

15 Since the difference between the boundary points of the selected group 

and the points outside and adjacent to the selected group (set to be zero) may 
be too large in value, the purpose of the time-domain window mask m[n] is 
to further modify the modified estimated CIR h[n] to reduce the difference 
of value between the boundary points of the selected group and the points 

20 outside and adjacent to the selected group. Thus, the determined initial value 
of the frequency domain TEQ impulse response will be more appropriate. In 
addition, in order to shorten the CIR h[n] by stepwise adapting the TEQ 108, 
the value of the "tail" (the last few points) of the generated target impulse 
response b[n] should be decreased gradually. In this manner, it is reasonable 

25 that the initial value of b[n] should have the similar property that the value 
of the "tail" of the initial b[n] should be decreased gradually. The purpose of 
the time-domain window mask m[n] is to further modify the modified 
estimated CIP h[n] not only to reduce the difference of value between the 
boundary points of the selected group and the points outside and adjacent to 

30 the selected group but to further modify the modified estimated CIP h[n] 
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such that the value of the "tail" of the initial b[n] is decreased gradually. 
Fig.4A and 4B show two examples of the time-domain window mask m[n], 
which can be mathematically represented by 



(Fig.4A) 



and 



0; others 



m[nl = 



n-k 



\k<n<i 



(i;others 



where i is an integer between k and k+ -1. 



It is emphasized here that the "shape" of the time-domain window 
mask is subject to modification according to system requirements and not 
limited to the examples of Fig.4A and 4B. 

10 In the step 34, h[n] and b[n] Bit transformed to the frequency domain 

by Fast Fourier Transform (FFT), which are represented by H(i) and B(i) 
respectively. 

B(i) determined in step 34 and W(i) determined in the steps 35 are the 
initial values of frequency-domain TEQ impulse response W[n] and target 
15 impulse response B[n], which can be applied in various adaptive TEQ 
algorithms. 

Fig.5 is a flow chart of the method for determining the initial values of 
frequency-domain TEQ impulse response W[n] and target impulse response 
B[n] according to the second embodiment of the present invention. The 
20 method includes steps of: 

50 producmg an estimated CIR h[n] according to the received 
symbols which are known by both of the receiver and the transmitter; 

51 calculating the total energy of each group of consecutive p -U 

points and selecting the group of consecutive v -Iw points with the 

s 



maximum energy, wherein Iw is the length of the TEQ impulse 
response; 

52 removing the group of consecutive v -Iw points with the maximum 

energy from h[n], combining the remaining N- v +lw points, and 

5 padding zero to the last v -K points of h[n] to produce a new N-point 

sequence h*[n]; 

53 transforming h*[n] into frequency domain to generate H'(i), 
i=O^N-l; 

54 determining W(i) which is the reciprocal of H'(i), i=0~N-l, 
10 wherein W(i) is the initial value of the frequency-domain TEQ 

iinpulse response; and 

55 transforming the original estimated CIR h[n] into frequency 
domain to generate H(i), and determining B(i) by multiplying W(i) 
with H(i), wherein i=0~N-l, and B(i) is the initial value of the 

15 frequency-domain target impulse response. 

Based on the method disclosed in this embodiment, the w[nj (i.e. the 
impulse response of the TEQ 108) is determined such that the convolation 
of w[n] and the new N-point sequence h'[n] generated in the step 52 is an 
ideal impulse response. In this manner, the valid length of the target impulse 
20 response b[n], which is generated by the convolution of the w[n] and the 
original h[n], will be shorter. Thus, the target impulse response b[n] 
determined by the method disclosed in this embodiment will be more 
appropriate that the valid length of the target impulse response will be equal 

to or shorter than that of the cyclic prefix, i.e. v points. 

25 In the step 51, the Locate maximum energy algorithm is performed. In 

the step 52, the group of consecutive v -1* points having the maximum 
energy are removed from h[n], the remaining points of h[n] are combined, 
and the last consecutive v -Iw points are set to be zero to generate a new 
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N-points sequence h'[n]. In the step 53, h'[n] is transformed into frequency 
domain by FFT, and H'(i) represents the FFT result of h'[n]. The W(i) and 
B(i) determined in the steps 54 and 55 respectively are the initial values of 
frequency-domain TEQ impulse response W[n] and target impulse response 
5 B[n], which can be applied in various adaptive TEQ algorithms. 

Fig.6 is a flow chart of a typical adaptive algorithm of the TEQ 108. As 
shown in Fig.6, the algorithm is performed in an iterative way to adapt the 
TEQ 108 to in order to shorten the CIR of the communication system. The 
algorithm includes steps of: 

10 60 determining initial value of frequency-domain TEQ impulse 

response Wk[i] and target impulse response Bk[i], wherein 
i=0-N-l ,and k is the adapting index representing k-th iterative 
operation; 

61 performing IFFT on Wic[i] and Bu[i] to generate Wk[i] and bi:[i] 
15 respectively. Then calculating the total energy of each group of 

consecutive points of Wk[i] and selecting the group of consecutive 
Iw points with the maximum energy, wherein 1^ is the length of the 
TEQ impulse response, and setting all the remaining points outside 
the selected group to zero to generate a modified time-domain TEQ 
20 impulse response W]£,w[i]. In the same time, calculating the total 

energy of each group of consecutive i/ points of bu[i] and 

selecting the group of consecutive i/ points with the maximum 

energy, and setting all the remaining points outside the selected 
group to zero to generate a modified time-domain target impulse 
25 response bk,w[i]- Then performing FFT on Wk,w[i] and bk,vv[i] to 

generate Wicw(i) and Bk,w(i) in frequency-domain, wherein i=0~N-l; 

62 performing Eic(i)= Bk,w(i)- WK,w(i) • H(i), i=0~N-l , wherein H(i) is 
the frequency-domain of h[n]; 

63 performing Wk+,(i)^ Wk,w(i)+ Ek(i) • H*(i), i=0~N-l, to adjust 
30 the value of the frequency-domain TEQ impulse response, wherein 
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H'(i) is the complex conjugate of H(i) and ii is a stepsize 
coefficient; 

64 performing IFFT on Wk+i[i] to generate Wb+i[i], then calculating 
the total ener^ of each group of consecutive Iw points of Wki.i[i] and 
5 selecting the group of consecutive U points with the maximum energy, 

wherein U is the length of the TEQ impulse response, and setting all 
the remaining points outside the selected group to zero to generate a 
modified time-domain TEQ impulse response Wk^i,w[i|- Then 
performing FFT on Wk+i.w[i] to generate Wt+i.wCi); 

10 65 performing Bw+i(i)=Wk*i,w(i) • H(i) to adjust the value of the 
frequency-domain target impulse response; and 

66 repeating the steps 61 to 65 until a predetermined converging 
condition is satisfied. 

In the step 60, it should be noted that the method to detennine the 
15 initial values of frequency-domain TEQ impulse response Wic[i] and target 
unpulse response Bk[i] is not limited to the methods disclosed in this > 
specification. 

In the step 62, Ek(i) can be regarded to be the error term generated by 
the adapting process. 

20 In the step 63, the equation Wicfi(i)= Wk.w(i)+ E)c(i) • H'(i) is the 

least-mean-square (LMS) algorithm, fi is the stepsize coefficient of the 

LMS algorithm. The value of /z may affect die performance of the LMS 

algorithm. If ii is too large, the LMS operation miay easily lead to 

divergence. If p. is too small, the LMS operation may lead to bad local 
25 maxima and the speed of convergence will be too slow. 

Based on the result of simulation, the value of fi is substantially in 
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proportion to the reciprocal of the power of H[i]. This condition may be 
shown by the following equation: fi = const / power(H[i]), wherein the 

value of const can be determined based on the result of simulation. It 
should be noted that divergence always occurs at an early stage of the 

5 adapting process if /z is too large. 

In a conventional approach, // is a time-invariant constant generated 
by experiments or simulations. However, in the embodiment of the present 
invention, /z is a time-variant fimction /z(k) throughout the adaptive 
process. Two methods for determining and controlling the value of the 
10 stepsize coefficient fimction p. (k) throughout the adapting process are 

disclosed in the specification. 

The first method is called the power-ratio method in the specification, 
the method can be represented by the following equation: 

const . power(fVk) 
^ power W pawer{Wk-Wk-\) 

1 5 Based on the above equation, // (k) will become large if power(Wk) is 

increasing and vise versa. The log function is used because the value of 
/?ow<rr(Ff*^^^^^^^^^^^^ increases exponentially during the adaptive 

process. If (Wv-Wt.,) is large, ^'"''''^^ypoy.erm-^-'^) ^^^^ ^^^^^ase 
such that divergence can be prevented. If (Wk-Wn) is small, 
2Q ^^^^^^^ypoyt,er(JVk -Wk -\) increase such that speed of convergence 

can be speeded up. According to the result of simulation, the value ofp (k) 

will increase gradually in the early stage of the adaptive process, and then 
approach to a constant in the late stage of the adaptive process. 

The second method is called the fitting curve method in the 
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specification. Through simulation, it is known that the value of /z (k) should 



keep small in the early stage of the adaptive process and then the value of /z 

(k) should become large in the late stage of the adaptive process. It can be 
accomplished by setting a predetermined fitting curve fimction. An example 

5 of stepsize coefficient fimction ii (k) is provided below, as shown in Fig.7: 

const 
poweriH) 

.wherein f(k)=k/M , when k<=M 

=1 , when k>M 

where k is the adapting index representing k-th iterative operation, and 
10 M is an integer which is preferably to be between 10 and 20. 

In the above example, ii (k) is kept small in the early stage of the 
adapting process (i.e. k<M), while in the late stage of the adaptive process 
(k>M), fi (k) is kept to be a large constant to prevent from converging too 

slow. It should be noted that the above example function is one of the 
15 embodiments and may not be the limitation of the present invention. 

In the step 65, the equation B,,-,i(i)=Wk*i,w(i) ♦ H(i) is called the 

zero-forcing criterion. The converging condition is not clearly determined in 
the step 66. Typically, the converging condition is predetermined to be that 
either the error Ek smaller than a predetermined threshold or perfonning the 
20 steps 61 - 65 for a predetermined period of time. 

While the present invention has been shown and described with 
reference to two preferred embodiments thereof, and in terms of the 
illustrative drawings, it should be not considered as limited thereby. Various 
possible modification, omission, and alterations could be conceived of by 
25 one skilled in the art to the form and the content of any particular 
embodiment, without departing from the scope and the spirit of the present 
invention* 
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